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Since the first analytical model of gravitational collapse proposed, using general relativity 
[1], many works have been done in this area. An important model is the one where the 
collapsing matter is represented by a massless scalar field and the space-time has a spherical 
symmetry (Einstein-scalar model). D. Christodoulou has pioneered analytical studies of 
that model [2]. In particular, he showed that for a self-similar scalar field collapse there 
are initial conditions which result in the formation of black holes, naked singularities and 
bouncing solutions. 

Many analytical results obtained by Christodoulou have been confirmed, by numerical 
studies of the Einstein-scalar model, by M. W. Choptuik [3]. Other authors have investigated 
the Einstein-scalar model numerically and analytically confirming the above results [4]. In 
particular, we may mention the works in references [5], [6]. There, the authors were able to 
reproduce qualitatively, in an analytical solution [7], the results derived numerically in [3]. 

The Roberts solution [7], has two arbitrary parameters a and (5. In the work in reference 
[5], the author fixed P — 1 and let a vary in the range [1/4, —oo]. Once /3 is fixed, a 
determines also the scalar field strength. It will be zero for a = 1/4 and it will increase 
for decreasing values of a. The scalar field eventually diverges for a — > — oo. Starting with 
a — 1/4, the scalar field is zero and one has Minkowski space-time. Then, for < a < 1/4, 
although the field is non- longer zero its strength is not enough to initiate the collapse. The 
space-time in this case represents a bouncing. For a = the gravitational collapse begins 
and a naked singularity is formed at r = 0. If one increases even further the field strength 
by letting a takes values in the range — oo < a < 0, one has the black hole formation. 
Therefore, if one turns the scalar field on and gradually increases its strength one will see 
all the collapse stages described above. Suppose now, one decides to turn the scalar field 
off, which will eventually happen for one has a limited amount of energy. One will see all 
the collapse stages described above, in the reverse order. Starting with the black hole and 
finishing with the Minkowski space-time. So, we may conclude that the Roberts solution 
does not produce an eternal black hole that will be left as a remnant of the scalar field 
collapse. Once the scalar field is turned off one loses all trace of the collapse. It is not 
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difficult to see that one may obtain the above result even if one let (3 vary in an appropriate 
range. 

One expects that the Schwarzschild black hole results from the collapse of a spherical 
symmetric neutral scalar field. So far, the only way to obtain this black hole as the result of 
the scalar field collapse was by joining an analytical solution to the Einstein-scalar equations 
that represents the collapse with the Vaidya solution [8]. Of course, this procedure generates 
the Schwarzschild black hole since it is a particular case of the Vaidya one. In reference [8] , 
the authors joined the Roberts solution [7] with the exterior Vaidya solution. 

In the present paper we would like to present a new solution to the Einstein's equation 
representing the self-similar, spherically symmetric, collapse of a minimally coupled, neutral, 
massless, scalar field, in four-dimensions. As we shall see this solution has two parameters, 
in the same way as Roberts' solution. Depending on the value of these parameters, we may 
have the formation of naked singularities and black holes. Since the black holes are identified 
as the Schwarzschild ones, we shall see how these black holes are produced as remnants of 
the scalar field collapse. 

We shall start by writing down the ansatz for the space-time metric. As we have men- 
tioned before, we would like to consider the spherically symmetric, self-similar, collapse of 
a massless scalar field in four- dimensions. Therefore, we shall write our metric ansatz as. 



where a{u, v) and r{u, v) are two arbitrary functions to be determined by the field equations, 
(ii, v) is a pair of null coordinates varying in the range (— oo, oo), and o?0^ is the line element 
of the unit sphere. 

The scalar field $ will be a function only of the two null coordinates and the expression 
for its stress-energy tensor T^p is given by [9], 



(1) 




(2) 



where , denotes partial differentiation. 
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Now, combining Eqs. (1) and (2) we may obtain the Einstein's equations which in the 
units of Ref. [9] and after re-scahng the scalar field, so that it absorbs the appropriate 
numerical factor, take the following form, 

2(rr,„^ +r,„r,J + e''^ = 0, (3) 



(4) 



(5) 



2(rV,„^ +rr,uv) = -r^{^,u^,v), (6) 

The equation of motion for the scalar field, in these coordinates, is 

r^,uv + r,u + r,v = . (7) 

The above system of non-linear, second-order, coupled, partial differential equations (3)- 
(7) can be solved if we impose that it is continuously self-similar. More precisely, the solution 
assumes the existence of an homothetic Killing vector of the form. 

Following Coley [10], equation (8) characterizes a self-similarity of the first kind. We can 
express the solution in terms of the variable z — v/u. 

In order to obtain our solution, we shall implement the self-similarity in the above equa- 
tions (3-7) in a slightly different way than previous works [5], [6]. We shall write the unknown 
functions as: ^{u,v) — ^{z), (t{u,v) — (t{z) and r{u,v) — r{z). It is clear that the main 
difference between our work and previous ones is in the way we write r. The system of equa- 
tions (3-7) will have the mentioned self-similarity if we re- write e^" in (1) in the following 
way. 
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where / = f{u, v) is an arbitrary function of the null coordinates u and v. f will also be 
written as f{z) when we implement the self-similarity. 

In terms of z and taking in account (9), the system of equations (3-7) becomes, 

2r(r' + zr") + 2z{r'y + zF^ = , (10) 

2zr" -%{F + 2zF') = - zr{^'f , (11) 
F 

2zr" + 2r' + ^[FF' + zFF" - ziF'f)] = - zri^'f , (12) 
F^ 

{zr^^y = 0, (13) 

where ' means differentiation with respect to z, f = F and equations (4) and (5) reduce to 
equation (11). 

The above system, equations (10-13), can be greatly simplified if we re- write it in terms 
of the new variable t = \nz. The resulting equations are, 

(r^)- = -e2*F2, (14) 
2r - 4r M + = -r6^ (15) 



2r + 2r (^J = -r6^ (16) 



(r2$)- = 0, (17) 

where the • means differentiation with respect to t. 

We solve the system of equations (14-17) by initially writing down a second order dif- 
ferential equation for r in the following way. First of all, we subtract equation (15) from 
equation (16) and manipulate it in order to find. 



where C is a real integration constant. Later, for a particular situation, we shall be able to 
identify the physical meaning of C. Then, we differentiate equation (14) once with respect 
to t and introduce, in the differentiated equation, the value of — e^*F^ from equation (14) 
and the information coming from equation (18). Finally, we integrate the resulting equation 
to find, 

r^r - 2Crf ^ B, (19) 

where S is a real integration constant. 

The physical meaning of B may be understood if we eliminate F and its derivatives from 
either one of the equations (15) or (16), with the aid of equation (18). They give the same 
result, which is. 

It is easy to see that equation (17) also leads to the same result above. 

Observing equation (20), we learn that 5 is a constant associated with the scalar field 
strength. More than that, we also learn that it cannot take positive values. Its domain is 
restricted to [0, — oo]. It is important to notice that equation (20) gives us a way to compute 
the scalar field once we know r as a function of t. 

In order to understand the meaning of the real constant C in equation (18), we shall 
consider the simphfied situation where B — 0. Prom equation (20) we see that $ is a 
constant. Without lost of generality we may set this constant to zero. It means that we are 
now dealing with a spherically symmetric vacuum solution to the Einstein's equations, so 
from Birkhoff 's theorem that solution is necessarily a piece of the Schwarzschild space-time 
[11]. 

We may identify the role played by C in the Schwarzschild space-time by first solving 
equation (19) for the present situation. This gives. 



r 2C 

t^- + —ln{Dr-2C) + E, (21) 
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where D and E are integration constants. Now, using the general expression of e^'^, which 
is derived from equation (9) and has the following value, 

e^^ = -^-^ , (22) 

uv 

we may compute g^^ for our present situation. With the aid of equations (14) and (21), we 

get, 

= m£L^. (23) 

uvr 

Finally, we may evaluate the coordinate independent quantity g'^/^r,ar,j3, for the present 
situation, and compare it with the same quantity evaluated for the Schwarzschild space- 
time. Using the values of r from equation (21) and g^^ easily derived from equation (23), 
we have, 

,^^r,.r„=i^L_J^. (24) 

Since the same quantity for the Schwarzschild space-time is 1 — 2M/r, where M is the 
Schwarzschild mass, we conclude that ii we set D — 1, C may be identified as the 
Schwarzschild mass. It is important to notice that this identification also tell us that C 
should not be negative in the general case {B,C ^ 0). If it were negative we would not 
be able to recover the Schwarzschild space-time as a limiting case of the general case. It is 
important to mention at this point that if one sets C — and B — one gets Minkowski 
space-time. 

With the aid of equations (21) and (23), for D = 1 and C — M, we may write the hne 
element equation (1) for the Schwarzschild solution, 

ds^ ^ - —(r - —) dudv + r\u, v)dn^ . (25) 



uv \ r 

One may transform this line element equation (25) to the usual one in terms of the (i, r) 
coordinates, 

= - (l - ^) + FTlM) + ■ P6) 
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In order to do that one has to use equation (21) and the coordinate transformation t — hi{uv). 

Let us consider the general case (S, M 7^ 0). Equation (19) cannot be solved analytically, 
as far as we know, therefore we solve it numerically. This equation can be written in the 
following suggestive form. 



Equation (27) reduces to the equation governing the case B — Q when we set r — > 00. This 
means that r{t) in the general case behaves asymptotically in the same manner as r{t) in the 
Schwarzschild space-time. Therefore, the space-times in the general case are asymptotically 
fiat. 

Then, taking in account this result we search for numerical solutions of equation (19) 
which asymptotically satisfy equation (21) and its derivative with respect to t, both with 
D — 1. Prom Eq. (21) we obtain that in the hmit r — > 00 we have r{t) — t and f{t) — 1. 
We have noticed that for a large number of different asymptotic values of r one obtains the 
same qualitative behavior of the solutions. In Figure 1, we exemplify this fact by tracing two 
graphics r x i for the space-time with M — 1 and B — —40, which we shall call 5'T140. In 
the first one. Figure 1(a), we set the asymptotic value of r to 10^. Therefore, the asymptotic 
conditions to solve Eq. (19) are r(lO^) = 10^ and r(lO^) = 1. In Figure 1(6), we set the 
asymptotic value of r to 10^. As a matter of definitiveness we have chosen in our calculations 
the same set of asymptotic conditions of Figure 1(a). 

For real, positive values of M and real, negative values of B we find an infinity set of 
solutions satisfying the above conditions. We shall call this set S. The r(i)'s, of the space- 
times in S, have the following behavior. They start collapsing from -|-oo initially linearly 
with t and then deviating from this regime until they reach r = 0, for certain values of t. 
These values of t depend in an unknown manner on B and M. In Figure 2(a), we can see 
r as a function of t for the space-time 5'T'140 and in Figure 2(6) for the space-time with 
M = 3 and 5 = -70, which we shall call ST370. 

From the curvature scalar R for the space-times in S, which is given by. 




(27) 
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D ^ (OR) 

^2(^2^2 + 2Mrr + B) ' ^ ' 

we can see that r = is a physical singularity. Performing a numerical investigation of a 
sufficiently large number of space-times in S we conclude that r = is the only physical 
singularity of these space-times. In Figure 3(a), we can see how R behaves as a function of 
t for the space-time 5'T140 and in Figure 3(6) for the space-time ST370 

Now, we may look for the presence of an apparent horizon in the solutions in S. We do 
that by investigating any change of sign in the quantity, g"^r,a This quantity has the 
below value for the space-times in 

'''"'■'^ " r2f 2 + 2Mrf + B ' ^^^^ 

Performing a numerical investigation of a sufficiently large number of space-times in S we 
conclude that they do not develop an apparent horizon. For large values of t, g"^r,a r,p 
equation (29) is a positive constant which has the value 1. As i decreases, it grows positive 
until it blows up at the r — singularity. Figure 4(a) shows the curve g°'^r,a T^p versus t for 
the space-time 5'T140 and Figure 4(6) for space-time ST'ilQ. We may also conclude, from 
the above result, that r = is a time-like singularity. 

Finally, we solve equation (20) for the space-times in S, together with the condition that 
$ asymptotically vanishes. Using the same asymptotic values of r and t, mentioned above, 
this condition translates to $(10^) = 0. Performing a numerical investigation of a sufficiently 
large number of space-times in S we find a scalar field with the following behavior. It grows 
negative from zero at t — > oo, initially very slowly, then more quickly as we diminish t. As 
we approach the physical singularity $ grows negative very quickly and finally blows up 
at the singularity. Figure 5(a) shows $ as a function of t for the space-time 5'T140 and 
Figure 5(6) for the space-time S'T370. It is not difficult to verify numerically that, in the 
asymptotic region, $ is and its derivatives with respect to t vanish as t — > oo. 

As a matter of completeness we mention that the particular case M = 0, S < 0, has 
an analytical solution with the same qualitative behavior of the numerical solution of the 
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general case (M > 0, B < 0). The main difference is that in the first case if one takes B — 
one obtains flat space-time instead of Schwarzschild space-time. Also, it is not difficult to 
see numerically that the so-called 'bouncing solutions', present in the collapse solutions of 
References [2], [3], [5], [6], can only appear from equation (19) for the unphysical situation 
of M > 0, S > 0. 

Collecting together all the above results we conclude that for the infinite number of 
space-times in S, parametrized by the constants M (M > 0) and B {B <0), ^ will collapse 
to r without the formation of an apparent horizon. Since r is a physical singularity 
for all elements of S, this collapse will give rise to naked singularities. These singularities 
are time- like and the space-times are asymptotically fiat. Combining together the above 
properties of these naked singularities with the fact that they define a region of non-zero 
measure in the parameter space of solutions, we can say that they qualify as possible counter- 
examples to the weak [12], [13] and strong [14] cosmic censorship conjectures. An issue that 
will be further investigated elsewhere. 

It is important to describe the space-times in S in terms of the original null coordinates 
{u,v). We start this program by noticing that, since t — ln{v/u) we have the following 
different ways to obtain the limit i — > oo if {u, v) are non-negative. They are: (i) — > oo 
while u remains finite and (ii) u while v is non-zero. We may also have the limit 
t — > oo if {u, v) are non-positive. Once they are physically equivalent to (i) and (ii) above, 
we shall restrict our attention to non-negatives {u, v) . Another important information is the 
location of r = in the (u, v) coordinates. As we have seen above, there is a certain t, let's 
say associated with the physical singularity. In terms of (xi, v), r = is the straight line 
passing through the origin with the equation: u — exp {—tojv. Collecting together the above 
information, we see that the space-times in S may be described as a piece of the m > and 
V >Q sector in (u, v) coordinates. Since, in the surface u — Q the space-times are fiat, we 
may extend them, smoothly, to the region it < as Minkowski space-times. Figure 6 shows 
the conformal diagram for a typical space-time in S. We may see that the space-time is 
divided, naturally, in two distinct regions. The first one is the Minkowskian region where 
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u < {!). Then, we have the collapse region where v > and u > (II). We can interpret 
this diagram in the following way: the scalar field starts collapsing from u = 0, the initial 
data surface in the present coordinates. The space-time is flat and $ and its derivatives are 
zero aX u — 0. Then, $ decreases with u until it blows up at the singularity r — 0. 

Finally, let us see one way to obtain the Schwarzschild black hole as the remnant of the 
scalar fleld collapse described by the solutions in S. Each solution in S depends on two 
real parameters M > and B <0. Prom equation (20) one can see that B determines the 
scalar field strength. Therefore, starting from Minkowski space-time M = B = Oatu — 
one turn the scalar field collapse on by increasing M and diminishing B. The scalar field 
strength increases and a time-like naked singularity is formed at r = 0. Suppose now one 
fixes the value of M = Mq and starts decreasing the scalar field strength by making B ^ 0. 
When B = the scalar field collapse ends and we are left with a Schwarzschild black hole 
of mass Mq. 

We would like to thank FAPEMIG for the invaluable financial support. 
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FIG. 3. Graphics oi R x t for space-times with different values of B and M. (a) ST140, (b) 
ST370. 
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FIG. 4. Graphics of g"'^r,a r,b x t for space-times with different values of B and M. (a) ST140, 
(b) ST370. 
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FIG. 5. Graphics oi ^ x t for space-times with different values of B and M. (a) ST140, (b) 
ST370. 




FIG. 6. Conformal diagram for a typical naked singularity solution. The null surface u = sep- 
arates the Minkowskian region (I) from the collapse region (II) where lies the time-like singularity 
at r = 0. 
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